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1 Introduction 



The method of Mellin-Barnes (MB) representation is one of the most powerful meth- 
ods to evaluate multiloop Feynman integrals. It is based on a very simple formula, 
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which is called MB representation and is applied to replace a sum of two terms raised 
to some power by their products in some powers. The contour of integration should 
be chosen so that the poles with a Y{. . . + z) dependence are to the left of the contour 
and the poles with a T(. . . — z) dependence are to the right of it. 

The simplest way to apply ^ is to represent massive propagators as continuous 
superpositions of massless ones, at the cost of introducing additional integrations over 
MB variables. More often, MB integrals are introduced at the level of Feynman or 
alpha parameters. Anyway, ([1]) is used in an appropriate way, with the goal to obtain 
integrals over loop momenta or Feynman/alpha parameters which can be taken in 
terms of gamma functions. As a result one obtains a multiple MB integral with an 
integrand involving gamma functions and powers of kinematic invariants. 

It is important to derive a MB representation for general values of the powers of the 
propagators. First, this provides unambiguous prescriptions for a choice of contours 
of integration. The resulting rule is similar to the rule for ([1]): when we understand 
a multiple MB integral iteratively and consider an integration over a variable 2, it is 
implied that the poles with a Y{. . . +Zi) dependence are to the left of the contour and 
the poles with a T(. . . — Zj) dependence are to the right of it. Second, such a general 
MB representation provides crucial checks and can be applied at any chosen values 
of the powers of the propagators. 

When evaluating a Feynman integral for specific powers of the propagators (in- 
dices), one starts from the general MB representation and obtains an integral of the 
form 
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where e = (4—d)/2 is the dimensional regularization parameter, aj, . . . , c^ - are rational 
numbers, are ratios of kinematic invariants and/or masses, and their exponents, 
dk, are linear combinations of e and z- variables. Typically, Cjj = ±1. 

Although results for dimensionally regularized Feynman integrals are practically 
needed in a Laurent expansion in e, one could try to evaluate a given MB integral 
(J2J) for general e. Then the e-expansion is performed in the result. However, for 
sufficiently complicated multiloop Feynman integrals, the evaluation at general e turns 
out to be impossible, and one proceeds with an e-expansion. 
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One of the advantages of this method is that the singularity structure in e can 
be resolved in a simple way. This procedure basically consists of taking residues and 
shifting contours, with the goal to obtain a sum of integrals where one can expand 
integrands in Laurent series in e. To do this one can apply two strategies formulated 
in pQ and [2] which we are going to call Strategy A and Strategy B, respectively. 

Strategy A is described and illustrated in numerous examples in Chapter 4 of 
[3]. It was applied, e.g., in [3J. According to strategy A, one performs an analysis 
of the integrand to reveal how poles in e arise. The guiding principle is that the 
product T(a + z)T(b — z), where a and b can depend on the rest of the integration 
variables, generates, due to the integration over z, the singularity of the type T(a + 
b). Indeed, if we shift an initial contour of integration over z across the point z = 
—a we obtain an integral over a new contour which is not singular at a + b = 0, 
while the corresponding residue involves an explicit factor T(a + b). This observation 
shows that any contour of one of the following integrations over the rest of the MB 
variables should be chosen according to this dependence, T(a + b). Hence one thinks 
of integrations in various orders and then identifies some 'key' gamma functions which 
are crucial for the generation of poles in e. Then one takes residues and shifts contours, 
starting from the first poles of these key gamma functions. The same analysis and 
procedure is applied to the contributions of the residues. 

Within Strategy B, one chooses an initial value of e and values of the real parts of 
the integration variables, Zi, . . . in such a way that the real parts of all the arguments 
of the gamma functions in the numerator are positive and one can integrate over 
straight lines. Then one tends e to zero and whenever the real part of the argument 
of some gamma function vanishes one crosses this pole and adds a corresponding 
residue which has one integration less and is treated as the initial integral within the 
same procedure. 

Strategy B is algorithmic in its character and, indeed, two algorithmic descriptions 
were formulated in [HE]. A public code called MB.m was presented in [6]. Strategy B 
was successfully applied, e.g., in [3 [8] and in many other papers. 

The purpose of this letter is to formulate Strategy A in a slightly modified form 
and to present the corresponding algorithm implemented in Mathematica. 

2 The modified Strategy A 

To present a modified Strategy A let us explicitly formulate what was implied in the 
initial Strategy A. When we take care of one of the key gamma functions we shift a 
contour and take a residue. Let T(Ai) with Ai = ai + hie + Ylij c ij z j be one of the 
key gamma functions in ([2]). Without loss of generality, we can consider e real. Then 
changing the nature of the first pole of this gamma function means changing the rule 
for an admissible contour, i.e. that, instead of the condition ReAj > when crossing 
the real axis in the process of the integration, we have the condition —1 < ReAj < 0. 
Let us denote this transition by replacing T(Ai) by T^(Ai). The initial rule for the 
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contour can be changed again and then we have the condition — n < Re A, < — n + 1 
for n — 2, 3, . . . with the notation r( n )(Aj). 

Within Strategy B one has straight contours in the beginning. Rather, in the 
modified Strategy A, we will be oriented at straight contours in the end. Apparently, 
it is desirable to achieve a minimal number of terms after the resolution of the singu- 
larities in e. To do this, let us try to search for contours which are going to have in 
the end of this procedure and for which the gamma functions in the numerator are 
changed, in the above sense, in a minimal way. 

To formalize this requirement, let us introduce the function cr(x) = [(1 — x) + ] 
where [. . .] is the integer part of a number and x + = x for x > and 0, otherwise. In 
other words, if — n < x < —n + 1 then u{x) = n for n > and cr(x) = for n < 0. 

So, let us set e = and search for contours, i.e. KeZi, for which the sum 



is minimal. 

After such a choice is done we identify gamma functions which should be changed, 
in the above sense, in order to arrive at a final integral where a Laurent expansion in 
e is possible. In fact, this step replaces the first step in the primary Strategy A where 
one identified such key gamma functions after the analysis characterized above. 

Then the second step in Strategy A is the same as in the old version: we take care 
of the distinguished gamma functions, i.e. take a residue and replace T by r^^Aj) 
(and, possibly, T^(Ai) by T^^Ai) etc.) We proceed iteratively, as in the previous 
strategy: every residue is considered from scratch, i.e. treated in the same way as the 
initial MB integral. 

Let us emphasize that although this strategy aims to minimize the number of 
resulting terms, we cannot exclude that there is another way of resolving the singu- 
larities in e that is the best one in this sense. For example, it can happen, in rather 
complicated examples, that different orders of changing the key gamma functions 
lead to different numbers of resulting terms. Still we believe that such a difference is 
negligible and that the Strategy A provides a resolution of the singularities in e at 
least very close to the theoretically best one. 

The difference of the new and the old Strategies A is minor. In the examples of 
[I] one can see that resulting contours were straight indeed. This difference can still 
be seen in the following simple example of the integral 



which, of course, can be evaluated at general e by the first Barnes lemma. Within 
the 'old' Strategy A, one can observe that there is no gluing of poles so that one 
can expand the integrand in e. However, a resulting contour cannot be chosen as 



Y ° ( Re ^4=o) = Y ° \ ai + Yl c ij Rez i 
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r(l + £ + z)r(-l/2 + e + z)r(3/2 - e - z)T(-z) dz 
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a straight line. Rather, within the 'new' Strategy A, we have to choose a gamma 
function to be modified and, as a result, we obtain a residue and an integral over a 
straight line which both can be expanded in e. 

3 Implementation in Mathematica 

The code should be loaded together with the package MB [6J by Czakon. It uses 
the MBresidue routine and some other private functions from that package; after 
obtaining the list of integrals one can continue the contour optimization and the 
numerical evaluation with the standard routines of MB. 

Formally the code can be described the following way: the initial function results 
in contour prescriptions (the arguments of T-functions have to be positive in the end). 
Now 

MBresolve[f unction, cont_pr] 

1. find a point with the minimal number of contour prescriptions broken. 

2. If this number is equal to zero, then return function as an answer 

3. / = list of broken contour prescriptions 

4. For each {x, n} in I do 

5. Try to evaluate tl = MBresolve[f unction, cont_pr'] 

where cont_pr' are obtained by replacing {x, n} with {x, n + 1} 

6. On exception move to the next For cycle 

7. represent x as ca + r, where c ^ 0, a is one of the integration variables 
and r does not depend on a 

8. Return tl - Sign[c]MBresolve[MBresidue[f unction, {a, —r/c}}, cont_pr"], 

where cont_pr" are inherited from cont_pr by removing {x, n} 
and performing the substitution a — > —r/c in all other functions 

9. End 

10. Throw exception 

The function should be a function of integration variables such that the poles are 
only due to the Gamma and PolyGamma factors (with linear arguments). The contour 
prescriptions cont_pr are a list of pairs x, n where a; is a linear function of integration 
variables and n is an non-negative integer. If n is equal to zero, such a term means 
that x has to be positive; if n is positive, then x has to be greater that — n and 
smaller than — n + 1. If the cont_pr parameter is missing, it is created automatically 
by considering all arguments of Gamma functions and their derivatives and assuming 
them all to be positive. An uncaught exception at the top level appears if there is a 
degenerate case and one is required to introduce extra regularization parameters. 

Hence, instead of using MBoptimizedRules [] or MBrules [] and, subsequently, 
MBcontinue [] , within MB.m [6j, one can now apply MBresolveG. 
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The corresponding Mathematica code MBresolve .m is public and can be found at 
http: //projects.hepforge.org/mbtools/| together with other tools for evaluating MB 
integrals. 

Here is an example of a tenfold MB representation derived loop by loop for the 
four-loop ladder massless on-shell diagram with pf = 0, i = 1,2,3,4, where s = 
(Pi + P2) 2 and * = (pi + Pzf- 

In[l] : 

« MB.m; 

« MBresolve. m; 

F = -(((-s)~(-5 - 4*ep - z7)*(-t)~z7*Gamma[l + zl]* 

Gamma [-1 - ep - zl - z2] *Gamma[-z2] *Gamma[-l - ep - zl - z3] * 
Gamma [-z3] *Gamma[l + zl + z2 + z3]*Gamma[2 + ep + zl + z2 + z3] * 
Gamma [zlO - z4] *Gamma[-zl + z4] *Gamma [-ep + zl + z2 - z4 - z5] * 
Gamma [-z5] *Gamma[-ep + zl + z3 - z4 - z6] *Gamma[-z6] * 
Gamma [1 + z4 + z5 + z6] * 

Gamma [1 + ep - zl - z2 - z3 + z4 + z5 + z6] *Gamma[-z7] * 

Gamma [1 + z7] *Gamma[-zlO + z7] *Gamma[-ep - zlO + z4 + z5 - z8] * 

Gamma [-z8] *Gamma[-ep + zlO - z7 + z8] * 

Gamma [-ep - zlO + z4 + z6 - z9] * 

Gamma [1 + ep - zlO + z7 - z8 - z9] *Gamma[-z9] * 

Gamma [-ep + zlO - z7 + z9]*Gamma[l + zlO + z8 + z9] * 

Gamma [1 + ep + zlO - z4 - z5 - z6 + z8 + z9] )/ (Gamma [-2*ep] * 

Gamma [1 - z2]*Gamma[l - z3]*Gamma[l - 2*ep + zl + z2 + z3] * 

Gamma [1 - z5]*Gamma[l - z6]*Gamma[l - 2*ep + z4 + z5 + z6] * 

Gamma [1 - z8]*Gamma[l - z9]*Gamma[l - 2*ep + zlO + z8 + z9])); 

Fcont = MBresolve [F, ep] ; 

Length [Fcont] 

MB 1.2 

by Michal Czakon 

improvements by Alexander Smirnov 

more info in |hep-ph/0511200| 

last modified 2 Jan 09 

MBresolve 1.0 

by Alexander Smirnov 

more info in arXiv: 0901. 03861 

last modified 4 Jan 09 

CREATING RESIDUES LIST 653.5156 seconds 

EVALUATING RESIDUES 15.2969 seconds 

Out [1] : 656 
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4 Discussion and perspectives 



Let us remind that, sometimes, a given MB representation can be ill-defined even for 
a well-defined Feynman integral, in the sense that a gluing of poles is present, i.e. 
one can distinguish a subproduct of the gamma functions such that the sum of there 
arguments, at general e, is equal to a non-positive integer number. To cure such a 
MB representation, one can introduce an auxiliary analytic regularization into some 
index, i.e. — > cii + y and then analytically continue the given MB integral, first, in y 
to the point y = and then e to the point s = 0. Let us stress that in such situation 
one starts, within the new Strategy A, with setting y = 0, e = and then searches for 
appropriate contours for which the gamma functions in the numerators are changed 
in the minimal way, similarly to the case without such analytic regularization. 

Let us emphasize that both Strategy A and Strategy B are based on the fact that 
singularities in e can be generated, because of gluing of poles of different nature, by 
the integration at compact regions. At least for planar diagrams, for which one can 
apply the code called AMBRE [9J based on the loop-by-loop strategy, this looks to be the 
only source of the singularities. However, for MB representations derived within the 
loop-by-loop strategy for nonplanar diagrams, there can be another source of polesjf] 
This feature can be exemplified by the massless two-loop nonplanar diagram with 
two external legs on-shelljj In the corresponding MB representation derived within 
the loop-by-loop strategy one meets, in particular, the following onefold MB integral 



There is no gluing of poles so that a pole in e cannot be generated by the integration 
over finite regions. Still a pole is generated and this can be seen by an explicit evalu- 
ation of this integral by closing the integration contour to the right and summing up 
the resulting series. This can be seen also by analyzing the asymptotic behaviour of 
the integrand at infinity. Setting z = x + iy and using formulae of the asymptotic be- 
haviour of gamma functions at large arguments in the complex plane, we can observe 
that the leading asymptotic behaviour when y — > +oo is l/y 1 ^ 26 which explains the 
appearance of the pole. Let us still mention that, for this concrete diagram, there is 
a better way to obtain a 'good' MB representation — see, e.g., Chapter 4 of [3]. 

For simple MB integrals, say, up to six- fold, the codes MB.m and MBresolve .m 
work more or less in the same way. For higher-fold MB integrals, e.g. ten-fold ones, 
MBresolve .m produces at first fewer integrals after the resolution of singularities in 
e. However, the command MBmerge of MB.m combines many terms, and the number 
of resulting integrals becomes smaller that produced by MBresolve .m. We should 
bear in mind, though, that in the current version of MB . m, this command combines 

3 See also a similar discussion in [8]. 

4 It was first calculated in an expansion in e up to e° in 10J . Now, results of expansion up to e 2 
[TT] . up to e 4 [12] and even a result for general e [T2] are available. 
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all integrals which have the same integration contours, while it is more natural to 
consider separately integrals with different patterns of gamma functions, both from 
analytical and numerical points of view. Presumably, the MBmerge command can be 
improved in this respect. 
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